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We present Shell Model Monte Carlo calculations for nuclei 
within the full major shell 50-82 for both protons and neu¬ 
trons. The interaction is determined solely by self-consistency 
and odd-even mass differences. The methods are illustrated 
for 124 Sn, 128 Te and 124 Xe. We calculate shape distributions, 
moments of inertia and pairing correlations as functions of 
temperature and angular velocity. Our calculations are the 
first microscopic evidence of y-softness of nuclei in this re¬ 
gion. 
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Nuclei with mass number 100 < A < 140 are believed 
"^to have large shape fluctuations in their ground states. 
Associated with this softness are spectra with an approx¬ 
imate 0(5) symmetry and bands with energy spacings 
ONntermediate between rotational and vibrational. In the 


geometrical model these, nuclei are described by potential 
energy surfaces with a minimum at f3 ^ 0 but indepen- 


8 ; 

f-A 

f-ylent of y [1], Some of these nuclei have been described in 
l/^terms of a quartic five-dimensional oscillator [2]. In the 
^(Interacting Boson Model (IBM) they are described by an 
(6) dynamical symmetry [3-5]. 

> Nuclei with 100 < ,4 < 140 fill the major shell between 
h 50 and 82 for both protons and neutrons, and conven¬ 
tional shell model calculations for the nuclei of interest 
3n the full space are impossible. However, with the intro¬ 
duction of Shell Model Monte Carlo techniques (SMMC j 
[6]. it has become possible to do exact calculations (up 
to a statistical error) in much larger model spaces [7] at 
zero and finite temperatures. This Letter presents the 
first fully microscopic calculations for soft nuclei with 
100 < A < 140 and compares them with the results of 
more phenomenological models. 

An important problem is the choice of the interaction. 
We used a pairing plus quadrupole interaction, where the 
pairing contains both monopole and quadrupole terms 
whose strengths are determined by odd-even mass dif¬ 
ferences. The quadrupole interaction is derived from a 
surface-peaked separable force [8] with a strength fixed 
by self-consistency. Such an interaction is expected to be 
a reasonable way of describing deformation and pairing 
phenomena. A major advantage of the interaction we 
construct is that it has a “good” Monte Carlo sign, so 
that calculations can be done quite accurately without 
resorting to extrapolation techniques [9], 


For the two-body interaction we take a multipole pair¬ 
ing force [10] plus an (isoscalar) self-consistent surface- 
peaked interaction. The latter is assumed to be separable, 
t>(r,r') = —\(dV/dr){dV/dr')S(f — r 1 ). where V(r) is the 
(spherical) mean-field potential. The angular delta func¬ 
tion is expanded in multipoles and only its quadrupole 
component is kept. We have 

Hi = ~ £ : • ft 

A /.< “ $ 


where :: denotes normal ordering and Pxy • Q p are pair 
and quadrupole operators given by 


II y II Jb)[a] a X ajjv 

ab 



In (2) a = ntj denotes a single particle orbit and dj m = 
( — ) J+m aj_ m . The quadrupole interaction strength is de¬ 
termined self-consistently. A change in the mean-field 
potential is related to a change in the one-body density 
p{ r) through 41 ir! = j dr'v(r,r')6p(r'). Using the in¬ 
variance of the one-bodv potential under a displacement 
of the nucleus, and the separable form of the two-body 
interaction, we obtain 


X 


-i 



2 dV dp 


dr r — -— . 

dr dr 


( 3 ) 


The spherical nuclear density in (3) is calculated from 
p(r) = (4tt)- ] J2afa R2 a( r )/ r< * sphere f a and R a are the 
occupation number and the radial wavefunction of orbit 
a. respectively, and the sum goes over both the core and 
the valence shells. The self-consistent result (3) gives 
\ = 0.053 MeV -1 fm 2 for 124 Xe. Since the interaction 
(1) is taken only in the valence shell, the value of \ found 
in (3) has to be renormalized. For a deformed harmonic 
oscillator single-particle potential, the core and valence 
quadrupole moments are equal [12], and the renormal¬ 
ization factor is 2. However, we use a Woods-Saxon po¬ 
tential V(r), and find that a renormalization factor of 
~ 3 is required to reproduce correctly the variation of 
the excitation energy of the first 2 + state in this regime. 

The single-particle energies are determined from a 
Woods-Saxon plus spin-orbit potential 
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hi = A + V(r) + A is0- s)V q — — + J z V c (r) . (4) 

2m r dr 2 

where V(r) = T o [1 + exp(r — i? 0 )/«] 1 with Bq = 
j’oA 1 ' 3 , and a Coulomb potential Vc of a uniformly 
charged sphere is included for a proton (r. = —1). We 
choose the paramctrization of Ref. [11], where T'o = 
—49.6 (l ± 0.86 N ^ Z ) MeV for protons and neutrons, re¬ 
spectively; a = 0.7 fin; /’o = 1.275 fin for protons and 
i d = 1.347 fm for neutrons. The spin-orbit potential 
has \( s = 0.488. j- 0 = 1.32 fin and the same diffuseness 
parameter as the central part. For 144 Xe, the resulting 
single-particle energies of 0 g 7/2 . ld 5 / 2 . Oh,, /2 , 2s*/ 2 . ld 3/2 
are -3.05. -3.35. -1.07. -1.04, -0.615 MeV for protons and 
-11.79, -12.08, -9.53, -10.21, -9.94 MeV for neutrons. No¬ 
tice that the Coulomb potential has the effect of placing 
the proton's h, ,/ 2 orbit below the d 3 / 2 and s,/ 2 orbits. 
When the central Woods Saxon potential is used for V(r) 
in ( 2 ). we find that the corresponding matrix elements of 
the quadrupole interaction in the proton single-particle 
basis differ by only a few percent from those in the neu¬ 
tron single-particle basis. We can therefore choose either 
of these sets. 

For the pairing interaction we include only monopole 
(A = 0) and quadrupole (A = 2) terms with g o = g 2 . 
go is determined in two steps. First, the pairing gap 
A is extracted from the experimental masses (or bind¬ 
ing energies B) of neighboring nuclei using [13] A„ = 

\ [B(N - 2, Z) - 3 B(N - 1. Z)+ 3 B(N.Z) - B(N + 1. Z)\. 
Next, using a particle-projected BCS approximation for 
the Hamiltonian (1), we find the value of g a that will re¬ 
produce the experimental gap for a spherical nucleus with 
the same mass number .4. For 144 Xe we find ga = 0.150 
MeV. The inclusion of quadrupole pairing is important 
in order to bring down the excitation energy of the 
state in the tin isotopes to about 1.3 MeV jwhich is 
2A ~ 2 MeV when only monopole pairing is included). 
Since both the monopole pairing and the quadrupole- 
quadrupole interaction are attractive, they satisfy the 
sign rule in the density decomposition and have a “good” 
Monte Carlo sign. The quadrupole pairing has compo¬ 
nents that violate that sign rule but they are all very 
small in comparison with the good sign components. Set¬ 
ting all the ‘‘bad" components to zero has no more than a 
591 effect on the spectrum. In identifying the “bad” com¬ 
ponents of the interaction, one should take into account 
that orbits of both parities are present in the 50-82 shell, 
and the modified sign rule should be used [9]. Thus, we 
are able to derive a quasi-realistic interaction that has a 
“good” Monte Carlo sign. 

We begin discussion of our results with the probabil¬ 
ity distribution of the quadrupole moment. This is an 
observable in principle, because it can be deduced from 
the expectation values of powers of the quadrupole op¬ 
erator. Observables in the SM MC are calculated by a 
weighted integration over their values for non-interacting 
nucleons moving in a fluctuating auxiliary field. We 
denote the single-particle evolution operator associated 


with an auxiliary field a by U a . and the quadrupole 
moment of each sample, (Q p ) a = Tr(Q fl U a .)/Tr(U a -), 
where Q p = r2 y 2p . The method used in a recent 
SMMC study of deformed nuclei [7] follows the prescrip¬ 
tion P{q,,) oc J2<t S ((Qi-‘)v ~ c hi)- Tllis method is valid 
provided Tr(Q^U„)/Tr(U c r ) = We can improve 

upon that prescription by explicitly forcing the quadratic 
moment to be correct. We calculate the variance of the 
Q operator for each sample. A 4 = Tr(Q 2 U,j)/Tr(U a ) — 
(Q)\ and assume that higher cumulants vanish. Then 


p (q ») lX J2 


5A ct 


exp 


HO!■)•■■ - r h ,) 2 


2A 2 


(• 3 ) 


In order for (5) to define a shape distribution we must 
associate a deformation a 2p with each quadrupole mo¬ 
ment {Q,,) a - Starting from a deformed nucleus, R = 
i?o(l + E„ ), we calculate its quadrupole moment 

by expanding its density p{ r) = p n (r — R) (where p 0 is 
the spherical density! to first order in deformation. We 
then find 


{Qij )<t — 4i?o 



a 2,i ■ 


( 6 ) 


po is calculated as below Eq. (3) but only within the va¬ 
lence shell. Eq. ( 6 ) is used to calculate from the 
measured We then transform a 2p , to the in¬ 

trinsic frame where a®o = 0 cos y and p 2 2 = ct 2 -2 = 
'3siny/\/2. The shape distribution P(3. 7 ) can be con¬ 
verted to a free energv surface through F(/3,^T) = 
— Flog [P(l 3,7)//J 3 I sin 3y |], where the unitary metric 
II„ /, = /3 4 | sin 3y | dSdgdBl has been assumed [14], 
The shape distributions of 14s Te and 144 Xe are shown 
in Fig. 1 at different temperatures. These nuclei are 
clearly 7 -soft, with energy minima at /3 ~ 0.06 and 
/3 ~ 0.15, respectively. Energy surfaces calculated with 
Strutinsky-BCS using deformed Woods-Saxon potential 
[15] also indicate 7 -softness with values of /3 0 compara¬ 
ble to the SMMC values. These calculations predict for 
124 Xe a prolate minimum with :3a ~ 0.20 which is lower 
than the spherical configuration by 1.7 MeV but is only 
0.3 MeV below the oblate saddle point, and for 14s Te a 
shallow oblate minimum with /3q ps 0.03. These 7 -soft 
surfaces are similar to those obtained in the 0 ( 6 ) sym¬ 
metry of the IBM, or more generally when the Hamilto¬ 
nian has mixed U( 5) and 0(6) symmetries but a com¬ 
mon 0(5) symmetry. In the Bohr Hamiltonian, an 0(5) 
symmetry occurs when the collective potential energy de¬ 
pends only on 3 [1]. Our results are consistent with a 
potential energy V([3) that has a quartic anharmonic- 
ity [ 2 ], but with a negative quadratic term that leads 
to a minimum at finite 3q. We have also estimated to¬ 
tal E 2 strengths from (Q 2 ) where Q = e p Q p + e„Q„ is 
the electric quadrupole operator with effective charges of 
e p = 1.5c and c„ = 0.5e, and extracted B(E 2;0 -4- 2j^) 
assuming that most of the strength is in the 2* state. We 
find B(E'2:Q -4 2j^) values of 663 ± 10, 2106 ± 15, and 
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5491 ± 36 e 2 fm 4 to be compared with the experimental 
values [16] of 1164. 3910. and 9103 e 2 fm 4 for 124 Sn. 128 Te 
and 124 Xe. respectively. Thus, the SMMC calculations 
reproduce the correct qualitative trend: the quantitative 
discrepancy would be ameliorated had we used a renor¬ 
malization factor similar to the potential field renormal¬ 
ization discussed earlier. 

Information on excited states can be obtained from 
strength functions. The energy centroid is given by 
E = S\ /Sq where S„ is the n-th moment of the relevant 
strength function, and can be calculated from the first 
logarithmic derivative of the imaginary time response 
function. This should be a reasonable estimate of the ex¬ 
citation energy of the collective state if the space of states 
does not include the giant resonance. We calculated the 
2+ excitation energy this way from the E2 response func¬ 
tion. The values of 1.12±0.02, 0.96A0.02 and 0.52±0.01 
MeV are in close agreement with the experimental val¬ 
ues of 1.2. 0.8 and 0.6 MeV for ]24 Sn. 128 To and 124 Xe. 
respectively. 

Another signature of softness is the response of the 
nucleus to rotations. We add a cranking field loJ z to 
the Hamiltonian and examine the moment of inertia as a 
function of the cranking frequency u>. For a soft nucleus 
we expect a behavior intermediate between a deformed 
nucleus, where the inertia is independent of the cranking 
frequency, and the harmonic oscillator, where the inertia 
becomes singular. This is confirmed in Fig. 2 which shows 
the moment of inertia I 2 = d{J z )lduj = ;3({J:) — {J.) 2 ) 
for 124 Xo and 128 Te as a function of oj, and indicates 
that 128 Te has a more harmonic character. The moment 
of inertia for lu = 0 in both nuclei is significantly lower 
than the rigid body value (« 43fi 2 /MeV for .4 = 124) 
due to pairing correlations. 

Also shown in Fig. 2 are { Q 2 ) where Q is the mass 
quadrupole, the BC'S-like pairing correlation (A'A) for 
the protons and (,7~). Notice that the increase in J> as 
a function of lu is strongly correlated with the rapid de¬ 
crease of pairing correlations and that the peaks in I 2 
are associated with the onset of a decrease in collectiv¬ 
ity (as measured by (Q 2 )). This suggests band crossing 
along the Yrast line associated with pair breaking and 
alignment of the quasi-particle spins at lu « 0.2 MeV 
((.J,) « 7 h) for 128 Te and w « 0.3 MeV ((J,) « 11 Ti) 
for 124 Xe. Our results are consistent with an experimen¬ 
tal evidence of band crossing in the Yrast sequence of 
124 Xe around spin of 10 h [17]. The alignment effect is 
clearly seen in the behavior of (J,) at the lower temper¬ 
ature which shows a rapid increase after an initial mod¬ 
erate change. Deformation and pairing decrease also as 
a function of temperature. 

We have analyzed the number of correlated pairs of 
these nuclei in their ground state. For a given angular 
momentum .7, we define the pair operators Aj M (ab) = 
I/yl — x a \,]' 1M ■ These operators are boson¬ 

like in the sense that they satisfy the expected com¬ 
mutation relations in the limit when the number of va¬ 


lence nucleons is small compared with the total num¬ 
ber of single-particle states in the shell. In the SMMC 
we compute the pair correlation matrix in the ground 
state ^2 M {A^j,u(ab)AjM{cd)), which is an hermitian and 
positive-definite matrix in the space of ordered orbital 
pairs ( ab ) (with a < b). This matrix can be diagonalized 
to find the eigenbosons B' nJM = Y. a b'*!■<*j( a b)A' JM (iah), 
where a labels the various bosons with the same angular 
momentum J. These eigenbosons approximately satisfy 

^( B jm§*) B $m( 0)) = n,j(a)S^ , (7) 

M 

where the positive eigenvalues n a j are the number of J- 
pairs of type a. We have calculated the total number of 
./-pairs (n,j = £ n a j) in the various pairing channels, 
and the results for the number of correlated pairs (after 
subtracting out the mean-field values [18]) are shown in 
Fig. 3. Since the number of neutrons in 124 Xe is above 
66 , they are treated as holes. For .7 = 0 and .7 = 2 we 
can compare the largest n n j with the number of s and d 
bosons obtained from the 0(6) limit of the IBM. In the 
latter we use the exact 77(6) formula [3] for the average 
number of pairs and multiply by the relative fraction of 
protons and neutrons to find the pair content for each 
type of nucleon. For 124 Xe the SMMC (IBM) results 
in the proton-proton pairing channel are 0.85 (1.22) s 
(.7 = 0) pairs, and 0.76 (0.78) d (.7 = 2) pairs, while in the 
neutron-neutron channel we find 1.76 (3.67) s pairs and 
2.14 (2.33) d pairs. For the protons the SMMC d to s pair 
ratio 0.89 is close to its 0(6) value of 0.64. However, the 
same ratio for the neutrons, 1.21, is intermediate between 
0(6) and SU( 3) (where its value is 1.64) and is consistent 
with the neutrons filling the middle of the shell. The total 
numbers of s and d pairs - 1.61 proton pairs and 3.8 
neutron (hole) pairs - are below the IBM values of 2 and 
6 , respectively. The fermion-based SMMC calculations 
indeed indicate pair correlations for higher .7 values. 

In conclusion, we have presented the first microscopic 
evidence of softness in nuclei in the 100-140 mass region, 
using the SMMC for the full 50-82 major shell. Work 
in progress on the systematic in this region would in¬ 
clude a Q p ■ Q ... force [5]. However, such calculations 
are more time-consuming as this interaction violates the 
Monte Carlo sign rule. 
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FIG. 1. Free energy surfaces for l2K Te (left) and l24 Xc 
(right) in the 8 — 7 plane at several temperatures. The con¬ 
tour lines are separated by 0.3 MeV and the lighter shades 
Correspond to lower energies. Notice the 7 -softness of the 
surfaces. 

FIG. 2. Observables for 124 X<g and 12 ‘ s Te as a function of 
cranking frequency mi and for two temperatures. J. is the 
moment of inertia. Q is the mass quadrupole moment, A is 
the J = 0 pairing operator, and J~ is the angular momentum 
along the cranking axis. 

FIG. 3. Number of correlated pairs of angular momentum 
J in 1 " 4 Xe. Different shades correspond to p-p, n-n, isoscalar 
and isovector p-n pairs. Neutrons are treated as holes. 
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